Astronomy & Astrophysics manuscript no. 1534 


February 2, 2008 


(DOI: will be inserted by hand later) 





X-ray source overdensities in Chandra distant cluster fields: a 
new probe to map the cosmic tapestry? 

N. Cappelluti\ M. Cappi\ M. Dadina\ G. Malaguti^ M. Branchesi^'^, V. D'Elia"^, and G.G.C. Palumbo^ 

^ lASF-CNR, Sezione di Bologna, Via Gobetti 101, 1-40129 Bologna, Italy 
e-mail: cappellutiObo . iasf . cnr . it 

Dipartimento di Astronomia Universita di Bologna, via Ranzani 1, 1-40127 Bologna, Italy 



2 



3 



IRA-CNR, Via Gobetti 101, 1-40129 Bologna, Italy 
* INAF-Osservatorio astronomico di Roma, Via di Frascati, 33 1-00040 Monteporzio Catone (Rome), Italy 
Received 06/25/2004 / Accepted 09/08/2004 

Abstract. The first systematic study of the serendipitous X-ray source density around 10 high z (0.24< z <1.2) 
clusters has been performed with Chandra. A factor ~ 2 overdensity has been found in 4 cluster fields, increasing 
to 11 the number of known source overdensities around high z clusters. The result is statistically highly significant 
(at >2(T per field, and >5a overall) only on scales 8'x8' and peculiar to cluster fields. The blank field fiuctuations, 
i.e. the most sensitive measurements to date of the X-ray cosmic variance on this scale, indicate la variations 
within ~ 15-25% around the average value. The first marginal evidence that these overdensities increase with 
cluster redshift is also presented. We speculate that the most likely explanation for what has been observed is 
that the x-ray sources are AGNs which trace filaments connected to the clusters. Nevertheless other possible 
explanations such as an "X-ray Butcher-Oemler" effect are not ruled out by the present results. If the association 
of the overdensities with large-scale structures and their positive relation with cluster redshift is confirmed, these 
studies could represent a new, direct, way of mapping the cosmic web and its parameters at high z. 

Key words. Galaxies: clusters:general - Galaxies: active - X-rays: galaxies: clusters - Cosmology: large-scale 
structure of Universe 



1. Introduction 

Recent studies have reported an unexpected excess of 
point-like X-ray sources on the outskirts of several dis- 
tant clusters of galaxies (Henry & Briel 1991, Cappi et al. 
2001, Pentericci et al. 2002, Molnar et al. 2002, Johnson et 
al. 2003). The source number densities measured around 
the clusters have been found to exceed the values observed 
in similarly deep observations of fields not centerd on clus- 
ters by a factor ^2. This is in contrast with optical and X- 
ray studies which suggest that AGNs are somewhat rare in 
cluster environments (Dressier et al. 1985, 1999; Branches! 
et al. in preparation) both in the nearby Universe and at 
high z. At least in one case (3C295 at z ~0.5), the excess 
number of sources has been demonstrated to be asym- 
metric with respect to the cluster, possibly indicating a 
filament of the large scale structure connected with the 
cluster (D'Elia et al. 2004). These pioneering studies em- 
phasize the potential of X-ray surveys for tracing the cos- 
Send offprint requests to: N. Cappelluti 



mic web at high z. The statistical significance and nature 
of these excesses however is not yet clearly established. 
The main limitations are the still poor statistics available 
from the X-ray studies, as well as the very few optical 
identifications available to date. In order to overcome the 
statistical problems a systematic and uniform analysis of 
the X-ray source population around 10 distant clusters 
of galaxies was conducted. The study was performed us- 
ing archival Chandra ACIS-I observations which are well 
suited to obtain X-ray source detections. Throughout this 
paper, we assume Ho= 70 km s^^ Mpc^^ and fl^l. Unless 
otherwise stated, errors are at the la level. 

2. Sample selection 

The cluster observations were selected from the Chandra 
archive using the following criteria: i) clusters should have 
z >0.2 so that a wide region around the cluster could be 
covered, ii) the cluster should be observed with the ACIS-I 
CCD (16' X 16') array with iii) the deepest exposure possi- 
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Table 1. Clusters and reference fields 



Field Name 


Obs. Id 


Exp. 
(ks) 


z 




Cluster fields: 










MS 1137+6625 


2228 


110 


0.78 


1.20 


CL J0848+4456 


927 


126 


0.58 


2.65 


RX J0910+5422 


2227 


124 


1.11 


2.09 


CL J2302+0844 


918 


110 


0.72 


4.91 


MS 2n,5.3-n44Q 


1667 


45 


0.58 


4.06 


Abell 2125 


2207 


81 


0.24 


2.74 


MS 1621+2640 


546 


30 


0.43 


3.57 


RDCS1252-29 


4198 


162 


1.20 


2.20 


CL J1113.1-2615 


915 


105 


0.72 


5.47 


IE 0657-56 


3184 


90 


0.29 


6.53 


Reference fields: 










CDFS 


582 


130 




0.80 


HDFN 


3389 


125 




1.60 


EXT-CDFS-1 


5015 


163 




0.80 


Groth-Wcst. Field 


4357 


85 




1.28 


Boot(-s-Finl(l 




115 




0.50 



*Galactic absorption column density, taken from Stark et al. 
1992, in units of xlO^" cm"^ 

a range of redshifts from 0.2 to 1.2 and with exposures 
from 30 ks up to 162 ks (see Table 1). 

These fields were compared with five blank fields (i.e. 
fields not centered on clusters, hereinafter called "refer- 
ence fields") of similar exposures (see Table 1), namely the 
"standard" and extended Chandra deep field south (CDFS 
and EXT-CDFS), the Hubble deep field north (HDFN), 
the Groth-Westphal strip and the "Bootes" field. In order 
to sample similar flux levels the analysis of the CDFS and 
of the HDFN was limited to observation Ids 582 and 3389, 
respectively. 

3. Data Reduction and Analysis 

Data reduction has been performed using the CIAO 
3.0.1 software package (http://cxc.harvard.edu) and 
applying the standard technique as follows: photon 
level 1 event lists were filtered so as to include 
only the standard event grades 0,2,4,6, the hot pix- 
els and columns were removed, the CTI (Charge 
Transfer Inefficiency) and gain corrections were applied 
(http: / /hca- www. harvard. cdu/'^alexey/acis/tgain), flick- 
ering pixels were removed and the intervals with back- 
ground rates larger than 3cr over the quiescent value were 
excluded. 

The analysis was performed in two energy bands 0.5-2 
keV (soft) and 2-10 keV (hard). The source detection was 
carried out chip-by-chip on the unbinned data running 
the wavelet algorithm in CIAO {wavdetect, Freeman et al. 
2002) on a scale of (V2)' pixel (i = 1, .., 8), i.e. 0.5" - -8". 
A threshold probability of spurious detections of 10~^ has 
been used (i.e. <1 spurious source per chip) and those 
sources with a detection significance <3a (see Johnson et 



choice the faintest detectable sources have at least ^ 8-9 
net counts, which corresponds to a minimum SNR of ^3. 
Following the indications given in Manners (2002), we esti- 
mated that this allows uncertainties due to the Eddington 
bias (Eddington 1940) to be lower than ~5%. The source 
net counts were corrected for the vignetting using an ex- 
posure map. In order to maximize the signal-to- noise of 
the sources detected in the hard energy band, the anal- 
ysis was first restricted to the 2-7 keV band and then 
count-rates were extrapolated to the 2-10 keV fluxes (e.g., 
Giacconi et al. 2001). The count-rate to flux conversion 
factor was obtained by assuming a power-law model with 
galactic absorption and photon index (F) estimated from 
the stacked spectrum of the detected sources in each fleld. 
The typical spectral index found is F ^ 1.3 — 1.4 in the 
soft band and F ^ 1.6 — 1.7 in the hard band. Detailed 
results on this procedure which gives a good estimate of 
the "true" average spectral shape of the detected sources 
is beyond the scope of this work and will be presented in 
a forthcoming paper. 

4. X-ray source number counts 

4.1. Sky coverage 

Due to instrumental effects, such as vignetting and PSF 
(Point Spread Function) off-axis degradation, the sensitiv- 
ity of the ACIS-I detectors varies significantly across the 
field of view. Therefore, to properly calculate the logN- 
logS of the fiux-limited samples under investigation, the 
sky coverage ft must be carefully estimated, ft defines the 
area of the detector, projected onto the sky, sensitive to 
a limiting fiux Sum and was calculated as follows. For 
each chip a background map was constructed from the 
image after excluding the detected sources and filling the 
remaining "holes" (using the dmfilth tool in CIAO) with 
the surrounding average background. The maps were then 
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Fig. 1. Sky coverage as function of flux in the MS 1137+6625 
fleld. The horizontal line represents the lower limit (20% of the 
total) of the sky coverage which gives the fiux limit used in the 
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Fig. 2. Reference fields logN-logS normalization vs. Slope relation in the soft {Left Panel) and hard (Right Panel) band. 
Circles indicate best-fit values obtained chip- by-chip (8'x8'), the thick diamond indicates the mean value error box calculated 
on scales of 8'x8' and cross indicates the mean calculated on 16'xl6' . 



rc-binncd into 64x 64 matrices of ~30" per pixel to smooth 
the local background variations and to sample the larger 
size of the ofF-axis PSF. Following Johnson et al. (2003), 
the background level in each bin was used to estimate the 
number of counts needed to obtain a >3cr detection (with 
wavdetect). The flux limit Sum (in units of erg cm~^ s~^) 
per image bin was then calculated using the formula: 

Sii^ = 3x B xV X f X {1 + ^/0.75 + {B x D^))/t (1) 

where B arc the background counts, V is the vignetting 
correction factor (given by the exposure map), / is the 
count-rate to flux conversion factor, D is the estimated 
source diameter and t is the exposure time. The sky cov- 
erage at a given flux Sum is then given by the sum of all the 
regions of the detector with a flux limit larger than 5iim 
and then converted to deg^. As shown in Fig. 1, the sky 
coverage rapidly decreases near the survey limiting flux 
for both the soft (continuous line) and hard (dashed line) 
energy bands. In order to prevent incompleteness effects 
only sources with flux greater than the flux corresponding 
to 20% of the sky coverage (see horizontal line in Fig. 1) 
were included in the analysis. 

4.2. logN-logS 

The cumulative source number counts N{> S) have been 
calculated chip-by-chip following the basic method de- 
scribed in Gioia et al. (1990) and using the following for- 
mula: 

n 

^(>^)=Ef|- deg-' (2) 

i—l 

where n is the number of detected sources and Cli is the 



A maximum likelihood method as given in Crawford 
et al. (1970) and Murdoch et al. (1973) was then used to 
flt the unbinned data with a single power law model of the 
form: 

Ni>S) = ki^r" deg-^ (3) 

where k and a are the normalization and slope of the 
power law. 

The only free parameter of the fit is a, while the nor- 
malization k is computed assuming that the predicted 
source counts at the flux limit are equal to what is ob- 
served. Errors in k were then calculated using the upper 
and lower limits of a. The value of k was calculated at 
flux levels of So=2xl0~^^ erg cm~^ s~^ in the soft band 
and So=lxlO^^^ erg cm^^ s^^ in the hard band; for the 
exposures considered here these correspond to the centers 
of the sampled flux intervals. 

5. Results 

The logN-logS of sources around the clusters has been cal- 
culated chip-by-chip, as well as over the whole (4 chips) 
FOV, and compared to the logN-logS of the reference 
flelds. 

5.1. Reference Fields 

As detailed above (see Sect. 2 and Table 1), 5 blank flelds 
(20 chips, i.e., a sky coverage of ~ 0.4 deg^) have been 
used as reference fields. Best-fit results of the logN-logS 
obtained chip-by-chip (8' x 8') are reported in Fig. 2. As 
shown in Fig. 2, the best-fit normalizations k are centered 
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Fig. 3. (Top) MS1137+6625 0.5-2 keV image, circles rep- 
resent the detected sources, the cross indicates the aim- 
point. {Bottom) 0.5-2 keV chip-by-chip logN-logS of the 
MS1137-I-6625 field. Dashed Une represents the 2cr confidence 
level of the CDFS-FHDFN logN-logS. 



5.2. Cluster fields 

As an example we first show the results obtained in the 
field of MS 1137-F6625. Fig. 3 (Top) clearly shows that 
the source density in the two boxed chips (one of which 
contains MS1137-I-6625) is about twice that in the other 
two chips. Chip-by-chip logN-logS have been calculated 
and are shown in Fig. 3. They are compared to the 
CDFS+HDFN logN-logS (CDFSN, ±2a deviations) that 
has a best fit set of parameters within 5% of the mean 
value given in Sect. 5.1. For two chips the logN-logS show 
a source surface density twice that in the CDFSN. 

As detailed above (see Sect. 2 and Table 1), a similar 
analysis was performed on 10 cluster fields for a total of 
40 chips, i.e. a sky coverage of ^ 0.8 deg^. Best-fit nor- 
malizations and slopes for the logN-logS in the soft and 
hard energy bands are shown in Fig. 4. and compared to 
the average values of the reference fields (marked with di- 
amonds in the figures). In six chip (from 4 cluster fields) 
excesses of a factor ~ 1.7-2 with a significance > 2o" are 
detected. These are clearly seen in Fig. 4 {Left Panels), 
and cluster- by-cluster statistics are given in Table 2. The 
excesses are located in tlic chip containing the cluster and 
in the adjoining chip. In the Lynx field (here named CL 
J0848-I-4456) the excess is observed in the only chip with- 
out a cluster (i.e. in this field three clusters were detected). 

The other 6 fields (24 chips) do not show any signifi- 
cant (> 2(7) excess. Nevertheless, in four out of these six 
cluster fields, given the low statistics, factor ~ 2 devia- 
tions cannot be statistically ruled out either. Altogether, 
we find at least 4 clusters out of 10 with significant over- 
densities on angular scales of 8' x 8' or 8' x 16'. 

It is worth noting that when calculated over the entire 
FOVs, these excesses are smoothed out and the logN-logS 
become statistically consistent with the reference fields on 
the same angular scales. 



Table 2. Statistical properties of the overdensities 



(standard deviation (7soft=70±3 deg^^) and khard=275tg^ 
deg~^ (crhard=36±5 deg~^), in the soft and hard energy 
bands, respectively. These average values (marked with di- 
amonds in Fig. 2) are in very good agreement with other 
recent studies (e.g. Moretti et al. 2003, Stern et al. 2002, 
Brandt et al. 2001, Giacconi et. al 2001, Rosati et al. 2002, 
Tozzi et al. 2001). 

Field to field variations, otherwise known as "cosmic 
variance" (see e.g. Soltan et al. 2001, Vikhlinin et al 1995), 
are within a spread of ~20-25% around these mean values. 
This is in agreement with the few other recent studies 
performed at similar fiuxes (e.g. Yang et al. 2003, Wang 
et al. 2004). Slope indeces are always consistent with the 
mean values of 0.63±0.13 in the soft and 0.75±0.15 in the 
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(1) 


(2) 


(3) 


(4) 


MS 1137-1-6625 


11 


2.47 


1.35% 


3.10 


0.20% 




I3t 


2.36 


1.83% 


2.57 


1.02% 


CLJ0848-h4456 


12 


2.25 


2.44% 


2.41 


1.60% 


RDCS0910-F5422 


I3t 


2.20 


2.78% 


2.97 


0.30% 




12 


2.03 


4.24% 


* 


* 


MS2053-0449 


10 


* 


* 


2.66 


0.78% 




13t 


2.23 


2.57% 


2.46 


1.43% 



(1) Significance of the excess in units of a (Soft). 

(2) Probability of upward Poisson fluctuation of reference field 

distribution (Soft). 

(3) Significance of the excess in units of a (Hard) 

(4) Probability of upward Poisson fluctuation of reference field 
distribution (Hard) . 

t Chip containing the cluster. 
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Fig. 4. logN-logS normalization vs. slope relations for overdcnse cluster fields in the soft (Top left panel) and hard band 
(Bottom left panel) . The same relation in cluster fields with non-significant or no excesses in the soft (Top right panel) and 
hard band [Bottom right panel). 




Fig. 5. Distribution of the logN-logS normalizations in cluster and reference fields, in the soft (Left Panel) and hard 
(Right Panel) energy bands. 
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Fig. 6. Overdensities (fcciusters/fcreference) as a function of clus- 
ter redshift, in the hard band. 



Fig. 5 shows the histograms of the normalization k 
in cluster and reference fields (chip- by-chip). The his- 
tograms clearly show that the distributions have an al- 
most Gaussian shape in the reference fields, both in the 
soft (Fig. 5, Left) and hard (Fig. 5, Right) energy bands. 
Instead, in the cluster fields, the distributions have a 
prominent tail at high values of k, containing about 20- 
25% of the fields. The significance of the tail can be quan- 
tified with a KS test which yields a probability that the 
two distributions are extracted from the same parent pop- 
ulation of 9.03x10"^ in the soft band and 2.55x10"^ in 
the hard band. 

5.3. Overdensities as a function of z 

A very interesting result obtained from the present anal- 
ysis is the apparent correlation found between the am- 
plitude of the overdensities and the cluster redshifts, as 
shown in Fig. 6. In this figure, the the (2-10 keV) k nor- 
malized to the value of the reference fields increases with 
the cluster redshift. The correlation is significant at a level 
of ~ 96% yielding a correlation parameter r=0.69 for 10 
independent points. The best-fit linear relation is here 
R=(0.63±0.24)z-F0.9±0.2 in the hard band while a simi- 
lar relation R=(0.55±0.10)z-M.06±0.16 is obtained in the 
soft band. 

6. Discussion 

The surface density of the X-ray serendipitous sources 
around 10 high z clusters has been investigated with 
Chandra and compared to the surface density of 5 ref- 
erence cluster-free fields observed with similar exposure 
times. Four of the cluster fields show a factor ^2 over- 
density of sources at a >2a significance level, while there 
clearly are no similar deviations present in the reference 



The fact that, as shown in Fig. 5, a fraction of ~25% 
of all fields near clusters shows overdensities demonstrates 
that this effect is unlikely to be associated to random sta- 
tistical fluctuations constitutes a peculiar property of clus- 
ter fields. With the present work, the number of known 
high z clusters showing such excesses increases to eleven^. 
At a first glance, this contradicts what is reported in Kim 
et. al (2004) who find no significant difference between 
cluster and cluster- free fields (62 fields in total). Those 
authors, however, searched for source count differences on 
larger angular scales^ than in the present work. This may 
have had the effect of smoothing any excess signal. In sup- 
port of the above argument we find that integrating our 
data over the whole (16'xl6') ACIS-I FOVs, we find no 
statistically significant difference between fields either. In 
other words, this is consistent with an increase, in clus- 
ter fields, of the angular correlation signal toward lower 
angular scales (see e.g. D'Elia et al. 2004). 

It should also be mentioned here that the present study 
yields one of the most sensitive measurements of the X-ray 
cosmic variance in empty fields available to date. Indeed, 
on scales of 8'x 8' as reported in Sect.5.1 and clearly shown 
in Fig. 5, we find (Icr) fiuctuations in number counts of 
the order of 17% and 13% in the soft and hard band, 
respectively. On scales 16'xl6', these values decrease to 
8% and 7% in the soft and hard band, respectively. 

Considering that, if the detected X-ray sources are as- 
sociated with the clusters themselves, their average lumi- 
nosity must be of the order of ~10^^~^'* erg/s at z=0.5 
and ~10^^~'''*'^ erg/s at 2;=1.2, they are most likely to be 
active galactic nuclei. This is also supported by the power- 
law shape of their summed spectra (see Sect. 3). Moreover, 
if connected with the clusters, the angular scales on which 
the overdensities emerge correspond to linear distances of 
the order of 3-7 Mpc (at z=0.5 and 1.2, respectively). This 
is somewhat greater than the typical dynamical radius of 
<3 Mpc for a massive cluster. 

A number of explanations have been proposed in the 
literature (e.g. Cappi et al. 2001) to explain this kind of 
excess of X-ray sources, e.g.: gravitational lensing, trigger- 
ing of AGN/starburst activity or filamentary clustering. 

The hypothesis that some of the detected sources are 
magnified by gravitational lensing seems rather unlikely 
because i) the shapes of the logN-logS do not get steeper 
at fainter fiuxes, as required by lensing in order to ob- 
tain a significant increase in the surface density of faint 
sources, and ii) the most distant clusters have Einstein 

^ The other 7 known cases being A2256 (Henry & Bricl 
1991), 3C295 and RXJ0030+2618 (Cappi ct al. 2001), 
MRC1138-262 (Pcntericci et al. ,2002), A2104 (Martini et al., 
2002) MS 1054-0321 (Johnson et al., 2003) and A1995 (Molnar 
et al. 2002). 

^ Indeed, Kim ot al. (2004) analyzed a total of 29 fields with 
clusters and 33 without and calculated the logN-logS either 
over a circular area per single FOV of radius 6.6 arcmin (i.e. 
with a sky coverage ~2 times larger than in this work), or 
adding together all fields reaching a cumulative sky coverage 
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angles of the order < 1', much smaller than the angular 
scale of the excess. This further confirms previous findings 
by Cappi et al. (2001) and Johnson et al. (2003), based 
on the calculations done by Refregier & Loeb (1997). 

Butcher & Oemlcr (1984) discovered that the fraction 
of blue (active) galaxies is very low (<<10%) in nearby 
clusters but increases linearly with the cluster redshift 
(percentages of several tens at z=0.5). This may well ex- 
plain, at least qualitatively, the redshift-dependence found 
here in the X-ray band (Fig. 6). Simulations suggest that 
this effect could be caused by the interaction of the intra- 
cluster with the interstellar medium in those galaxies 
which are falling into the cluster potential well (see e.g 
Evrard et al. (1991)), thereby triggering AGN/Starburst 
activity. However, such an explanation faces the problem 
that the scales involved here (>3 Mpc) are substantially 
greater than typical virial radii, thus it is not clear how 
much diffuse gas would be available at those distances to 
trigger such an interaction. 

The third hypothesis, that the excesses could repre- 
sent filaments of the large-scale structure of the Universe, 
is in our opinion the most likely explanation. In the cur- 
rent vision of the Universe, clusters lie in the knots of a 
filamentary web-like network of galaxies and Dark Matter 
(see e.g. Peebles, 1993). This has been demonstrated by 
both hydrodynamic simulations (Bond et al. 1996,Colberg 
et al. 2000) and by extensive optical imaging and spectro- 
scopic observations of galaxies around clusters (e.g. Gioia 
et al. 1999, Ebeling et al. 2004, Le Fevre et al. 2004). 

If AGNs trace galaxies (as observations clearly indi- 
cate, see for example Fig. 13 in Kauffmann et al, 2004), 
it is thus possible that we are observing AGNs in fila- 
ments connected with the clusters themselves. The de- 
tailed study by D'Elia et al. (2004) of the field surround- 
ing the cluster 3C295 at ^=0.5 clearly shows a strong 
and asymmetric clustering of X-ray sources on scales of 
a few arcmin. This dataset is, to our opinion, one of the 
best pieces of evidence to date for a filament at z ~0.5 
traced by X-ray emitting AGNs. Spectroscopic identifica- 
tion of the sources is being carried out and some redshifts 
close to that of the cluster give a first hint of a spatial 
connection between the overdensity and the cluster itself 
(D'Elia, private communication). Ohta et al. (2003) re- 
cently discovered in the optical band an overdensity of 
bright QSOs (a QSO-supercluster) and several groups of 
galaxies in the cluster field CL J0848-I-4456 at z~ 1.2. 
This cluster field also shows an excess of X-ray sources 
(Fig. 4), and thus could be related to the large-scale struc- 
ture at z ~1.2. Unfortunately, these optically identified 
quasars are spread over an area (30'x30') much larger 
than the Chandra FOV (see Fig. 3 in Ohta et al. 2003). 
Therefore only four identifications (one QSO at z ~1.26, 
one at z ~0.88 and two at z ~0.58) are in common mak- 
ing conclusions uncertain. Assuming that the excesses are 
indeed linked to the large scale structure of the Universe, 
the redshift-dependence of the overdensities (Fig. 6) could, 



In order to firmly establish whether the X-ray overden- 
sities found here are associated to large-scale structures 
at high z, extensive spectroscopic follow-up observations^ 
similar to those performed by e.g. Gioia et al. (1999), 
Ebeling et al. (2004), are required. 

But clearly X-ray studies of source overdensities 
around high z clusters might be a rather efficient way of 
mapping the cosmic web at high z. 

7. Conclusions 

The first systematic study of the serendipitous X-ray 
source density around high z galaxy clusters has been 
performed with Chandra. For the first time the results 
of study have been compared to a sample of cluster-free 
fields which also allowed a detailed measurement of the 
cosmic variance. Forty percent of the cluster fields present 
strong statistically significant overdensities. The first, al- 
though not robust, evidence for a redshift dependence of 
these excesses has also been presented. 

We speculate that filaments connected to the clus- 
ters are the most likely astrophysical interpretation of the 
above results. However, thought less likely, other effects 
such as an "X-ray Butcher Oemler effect" cannot be com- 
pletely ruled out. 
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